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ABSTRACT 

Stellar models with masses ranging from 0.5 to 1.3 Mq were constructed in order to 
compare to young cluster observations of Li and of rotation velocities. The amount of Li 
depletion in cool stars is sensitive to the amount of overshoot at the base of the surface 
convection zone, and the exact metallicity of the models. Even when this uncertainty is 
taken into account, the Li observations are a severe constraint for the models and rule 
out standard models, and pure diffusion models. Stellar models which include diffusion 
and rotational mixing in the radiative regions of stars are able to simultaneously match 
the Li abundances observed in the Pleiades, UMaG, Hyades, Praesepe, NGC 752 and 
M67. They also match the observed rotation periods in the Hyades. However, these 
models are unable to simultaneously explain the presence of the rapidly rotating late G 
and K stars in the Pleiades and the absence of rapidly rotating late F and early G stars. 
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1. Introduction 

Standard stellar evolution models have been very 
successful in explaining the morphology of the Hertz- 
sprung-Russell diagram. However, as new types of 
observations have become available, it has become 
clear that the standard models are missing a key in- 
gredient - the fact that stars transport material in 
stellar radiative regions. The evidence for particle 
transport in the radiative regions of stars comes from 
observations of ^Li in open cluster stars^ (see Pin- 
sonneault 1994). ^Li is an excellent tracer of particle 
transport in stellar radiative regions as it is easily 
destroyed at temperatures of ~ 2.6 x lO^K. Thus, 
there exists a radius below which no ^Li exists. If 
this radius lies below the surface convection zone of 
a star then the ^Li observed at the surface of the 
star should not change with time. Recent observa- 
tional studies have provided convincing evidence that 
main sequence ^Li depletion occurs in moderate mass 
(M ^ 1.3 Mq) stars (Hobbs & Pilachowski 1988; 
Boesgaard 1991; Thorburn ei al. 1993; Soderblom et 
al. 1993a). Significantly, beryllium is also observed to 
be depleted in the Sun (Anders and Grevesse 1989 and 
references therin) and in some field stars (Boesgaard 
1976). Because beryllium burns at higher tempera- 
tures than ^Li , this indicates that particle transport 
takes place throughout the envelope (albeit at a re- 
duced rate at increased depth). 

Main sequence ^Li depletion is in clear contradic- 
tion with standard stellar evolution models. Several 
mechanisms have been advanced to explain this deple- 
tion: main sequence mass loss, microscopic diffusion 
of ^Li out of the surface convection zone, and mix- 
ing driven by rotation. Swenson & Faulkner (1992) 
found that ^Li depletion induced by main sequence 
mass loss has serious problems when applied to clus- 
ter observations. Thus, the ^Li observations provide 
compelling evidence that some form of particle trans- 
port occurs in the radiative regions of stars. 

Evidence that rotation may be important in deplet- 
ing ^Li comes from detailed studies of the Pleiades 
(Soderblom et al. 1993a), the Hyades (Thorburn et 
al. 1993) and the Praesepes (Soderblom et al. 1993b) 
which found that stars with the same mass, compo- 
sition and age could have different ^Li abundances. 
Both standard stellar evolution models, and models 



"^In stars, Li may exist in two different isotopes: ^Li and ^Li. ft 
is extremefy difficuft observationaUy to distinguisfi between ^Li 
and ^Li, tfius tfie observations typicaffy determine tfie totaf Li 
present in a star, fn meteorites, 7.5% of tfie Li atoms are ^Li 
(Anders & Grevesse 1989). Tfius, it appears tfiat tfie totaf Li 
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Li observations actuaffy measure ^Li. 



which include diffusion predict that the Li depletion 
in a star is a unique function of age, mass and metal- 
licity. Some other property must be important in de- 
termining the amount of particle transport within the 
radiative region of a star. A logical candidate for this 
property is the initial rotation velocity, for it is ob- 
served that T-Tauri stars with similar masses have a 
wide range in rotation velocities (Bouvier 1991). The 
fact that rotation can induce mass motions in stel- 
lar radiative regions has been known for some time 
(Eddington 1925; Sweet 1950; see Zahn 1993 for a 
recent review), however these motions have typically 
been ignored in stellar evolution calculations (but see 
Charbonneau & Michaud 1991; Charbonneau 1992). 

Pinsonneault et al. (1989) and Pinsonneault, Ka- 
waler & Demarque (1990) presented detailed stellar 
evolution models which included the effects of rota- 
tional mixing. Our philosophy is to treat the prob- 
lem of stellar models with rotation as an initial value 
problem, in the spirit of Endal and Sofia (1976; 1978; 
1981). We specify the initial conditions, and a set of 
physical assumptions about angular momentum loss, 
transport, and particle transport. We then calibrate 
our models on the Sun, requiring that the surface ro- 
tation rate and lithium depletion match the solar val- 
ues. In practice, this determines the constant in the 
angular momentum loss law and the ratio between the 
diffusion coefficients for angular momentum transport 
and mixing. We then apply the same set of physics to 
models with a range in initial mass and angular mo- 
mentum and explore the observational consequences 
of changes in the physical assumptions. For simplic- 
ity, we have assumed that parameters in our models 
are the same (and constant) regardless of the mass or 
age of the star being modeled. The situation inside 
an actual star may be considerably more complicated. 

The surface rotation rates and light element abun- 
dances of low mass stars as a function of mass and 
age are sensitive to changes in the input physics and 
initial conditions. The observations allow us to rule 
out some classes of models, and can also be used as 
a guide to which uncertainties in the models have a 
significant impact on observable quantities. In pre- 
vious work, we have studied rotational models of the 
Sun (Pinsonneault et al. 1989) and open cluster stars 
(Pinsonneault et al. 1990). We begin by summarizing 
previous results, and then discuss the consequences of 
including microscopic diffusion and our improvements 
in the input physics. 

Angular momentum loss drives the rotational evo- 
lution of low mass stars. A general expression for 
the angular momentum loss rate is a constant multi- 
plied by the surface angular velocity to some power. 
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The exponent affects the time dependence of the spin 
down at late ages (Kawaler 1988), while the constant 
is adjusted to reproduce the Sun. Efficient angu- 
lar momentum loss in rapid rotators efficiently spun 
them down during the pre-MS, resulting in maximum 
ZAMS rotation rates of order 50 km/s at 30 Myr 
and 20 km/s at 70 Myr for solar mass models; this 
is in contradiction with observations in young stars 
(Stauffer and Hartmann 1987). A wide range in ini- 
tial angular momentum was reduced to a small range 
in surface rotation in ~ 2 x 10® yr, in agreement with 
rotation periods in the Hyades (Radick et al. 1987). 
These properties persisted for a wide range in param- 
eter space. 

Models with different initial angular momentum 
experience different degrees of rotational mixing; the 
time scale of angular momentum transport affects the 
timing of the mixing but not the overall amount in so- 
lar age systems because of the solar calibration. The 
timing of the lithium depletion was found to be con- 
sistent with the open cluster data in solar analogues, 
but the models of late F stars in the Hyades were 
over-depleted in Li with respect to the data - indicat- 
ing potential problems with the mass dependence of 
the particle transport. The small data sets available 
at the time, did not allow for a quantitative com- 
parison of the predicted and observed dispersion in 
lithium caused by mixing, although the existence of 
a dispersion could already be seen in some systems. 

In these models which included rotational mixing, 
but not diffusion, the sensitivity of transport to gra- 
dients in the mean molecular weight were found to 
have only a small impact on the surface properties 
because in standard stellar models such gradients are 
only large in the deep core, which is far from the re- 
gion where lithium is destroyed and has only a small 
moment of inertia. Since the work of Pinsonneault et 
al. (1989; 1990) was completed it has become clear 
that microscopic diffusion could have an important 
impact on stellar models (e.g. Proffitt & Michaud 
1991). In addition, relatively easy to use formulae 
for the microscopic diffusion coefficients have become 
available (Bahcall & Loeb 1990; Michaud & Proffitt 
1993). For this reason, a study of the combined ef- 
fects of rotation and diffusion was undertaken. In 
addition, the opportunity was taken to update the 
physics (both rotational and standard) used within 
the code. In particular, a saturation level to the an- 
gular momentum loss law has been added, in order to 
allow for a wide range of rotation velocities near the 
ZAMS. These changes have had a significant impact 
on the models. 

Microscopic diffusion leads to the separation of hy- 



drogen and helium, naturally leading to sharp gradi- 
ents in mean molecular weight at the base of the con- 
vection zone (see Vauclair 1983; Michaud and Char- 
bonneau 1991 for good discussions of the underlying 
physics). Microscopic diffusion also causes ^Li to sink 
with respect to hydrogen in stellar envelopes (see Vau- 
clair 1983, Michaud and Charbonneau 1991 for good 
discussions of the underlying physics) when it is fully 
ionized. For surface convection zones which are suf- 
ficiently thin, ^Li can be partially ionized and radia- 
tion pressure can cause ^Li to rise. Michaud (1986) 
proposed diffusion processes as the explanation for 
the strong depletion of ^Li seen in mid-F stars in 
the Hyades by Boesgaard and Trippico (1986). The 
timescale of diffusion increases with increased surface 
convection zone depth, making diffusion an unlikely 
explanation for later-type stars. 

Chaboyer, Demarque & Pinsonneault (1994, Pa- 
per I) examined the effect of varying the sensitiv- 
ity of rotational transport to gradients in the mean 
molecular weight. This degree of sensitivity was pa- 
rameterized by , large values of imply that the 
rotational mixing is strongly inhibited by gradients 
in mean molecular weight. In Paper I we found that 
fp. ^ 0-1 is required for appreciable main sequence 
^Li depletion to occur. This is due to the fact that 
diffusion can lead to a sharp gradient in mean molecu- 
lar weight at the base of the surface convection zone, 
which effectively halts the rotational mixing to the 
surface. 

In this paper, we will examine the ^Li depletion 
due to diffusion and rotational mixing and compare 
to observations of ^Li abundances and surface rota- 
tion velocities in young cluster stars. Charbonneau 
& Michaud (1991) and Charbonneau (1992) studied 
the combined effect of diffusion and advection due 
to meridonal circulation (Charbonneau & Michaud 
1991; Charbonneau 1992). Charbonneau & Michaud 
treated the meridonal circulation as a advection pro- 
cess; we treat meridonal circulation as a diffusion 
process. This is in agreement with the work of 
Chaboyer & Zahn (1992), who showed that appro- 
priate rescaling of the turbulent diffusion coefficient 
allows the combined advection/diffusion process in- 
duced by meridonal circulation to be treated solely 
as a diffusion process. This work includes the effects 
of angular momentum loss, as well as many rotational 
instabilities which was not studied by Charbonneau 
& Michaud. 

This paper is a follow up to Paper I which exam- 
ined the interaction of diffusion and rotation in the 
Sun. Observations of young cluster stars provide an 
excellent test of stellar models. Stars in a given cluster 



3 



have similar ages, distances and metallicities. These 
quantities are known far more precisely for cluster 
stars than for field stars. By studying a single clus- 
ter, it is possible to obtain information regarding the 
mass dependence of the ^Li depletion and surface ro- 
tation velocities. By comparing clusters of different 
ages, it is possible to determine how the ^Li depletion 
and surface rotation velocities vary with time. For 
this reason, a comparison between models and obser- 
vations in young clusters provides a crucial diagnostic 
of particle transport which occurs in the radiative re- 
gions of stars. Our angular momentum loss law is in- 
applicable to early-type stars, and we do not include 
radiative levitation in our models; for these reasons, 
we concentrate on models with masses of 1.3 Mq and 
lower. Hence, we will concentrate on comparing our 
models to the observations of ^Li abundances and ro- 
tation velocities in young stars with M < 1.3 Mq. 
A related paper, Chaboyer & Demarque (1994) ex- 
amines the constraints that the observed halo star 
^Li abundances put on particle transport in stellar 
radiative regions. 

Although cluster observations contain a great deal 
of information regarding the timescale and mass de- 
pendence of ^Li depletion, the Sun provides a unique 
calibration of stellar models since its properties (age, 
chemical composition, luminosity) are known far more 
accurately than any other star. This allows one to cal- 
ibrate the mixing length parameter which is used to 
model convection. By comparing photospheric and 
meteoritic abundances it is possible to determine the 
amount of destruction of ^Li. The Sun is the only 
star for which it is possible to make this measurement. 
For this reason, the Sun offers a unique opportunity 
to calibrate some of the parameters which are present 
in stellar models. Hence, all of the models we present 
in this paper will be calibrated to match the observed 
solar parameters, as discussed in Paper L Full details 
of the stellar model construction, including descrip- 
tions of the particle and angular momentum trans- 
port and angular momentum loss law are found in 
Paper L The mixing parameters and the properties 
of our calibrated solar models are shown in Table 1 
of Paper L Each set of mixing parameters was given 
a two letter code which will be used to refer to them. 
In these calibrated models, the mixing length, solar 
helium content and total metallicity have been var- 
ied in order to match the observed solar luminosity, 
radius and surface heavy element to hydrogen ratio. 
In addition, models which include rotation have been 
calibrated to yield the present solar surface rotation 
rate and ^Li depletion factor by varying the efficiency 
of mixing of chemical species, and an multiplicative 



constant in the angular momentum loss law. 

In §2. we provide a brief description of the obser- 
vational database which is used to test the models. 
The construction of stellar models and isochrones is 
discussed in §3. where some standard stellar mod- 
els are presented and shown to disagree with the 
^Li observations. Models which only include micro- 
scopic diffusion are presented in §4. In §5. the com- 
bined rotation and diffusion models are presented and 
found to be in reasonable agreement with the young 
cluster observations. Finally, the main results of this 
paper are summarized in §6. 

2. The Observational Database 

Observations of ^Li and surface rotation velocities 
in young clusters have the potential to delineate the 
time and mass dependence of angular momentum and 
mass transport within stars. Due to the advent of 
high efficiency CCD spectrographs, the database of 
cluster stars with observed ^Li abundances and ro- 
tational velocities has grown enormously in the last 
few years. The Pleiades, Hyades and Praesepe clus- 
ters have been extensively studied. The Pleiades 
has an age of 70 Myr which is particularly inter- 
esting as it implies that the solar mass stars have 
just arrived on the main sequence. Hence, obser- 
vations in the Pleiades give us information on the 
pre-main sequence ^Li depletion and angular momen- 
tum loss. The Hyades and Praesepe are ten times as 
old, 700 Myr (Barry, Cromwell, & Hege 1987) and 
so yield information on processes which occur on the 
main sequence. One complicating factor is that the 
Pleiades has a metallicity of [Fe/H] = -0.034 ±0.024 
(Boesgaard & Friel 1990), while the Hyades is more 
metal rich, [Fe/H] = -^0.127 ± 0.022 (Boesgaard & 
Friel 1990). Praesepe has a metallicity of [Fe/H] = 
+OmS ± 0.039 (Friel & Boesgaard 1992). In the 
Pleiades, high quality ^Li observations exist for 128 
stars (Soderblom ei al. 1993a). Accurate t; sin i de- 
tections have been made for 89 stars and upper lim- 
its of 7 — 10 km/s exist for 61 stars (Soderblom et 
al. 1993c; we have excluded stars with high redden- 
nings). In the Hyades, a uniform set of high quality 
^Li abundances for 52 stars and upper limits for 13 
stars has been obtained by Thorburn et al. (1993). 
The Hyades are unique among the clusters in that ro- 
tational periods exist for 22 stars (Radick et al. 1987). 
The fact that we do not have to worry about inclina- 
tion effects makes the use of rotation periods far supe- 
rior to t; sin i observations. Accurate ^Li abundances 
have been obtained for 63 stars in the Praesepe by 
Soderblom et al. (1993b). The Praesepe cluster has 
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an age very similar to the Hyades, so it is interesting 
to inter-compare these two clusters. 

In addition to the large data bases for the Pleiades, 
Hyades and Praesepe, there are a few other clus- 
ters which have interesting ages for which some ob- 
servations are available. There are no nearby clus- 
ters with an age of 200 to 300 Myr; however we 
are in a midst of a group of stars whose common 
origin can be discerned from their common kine- 
matics and chromospheric emission. The Ursa Ma- 
jor Group (UMaG) has an age of 300 Myr (Tomkin 
& Popper 1986, Soderblom & Mayor 1993) and a 
metallicity of [Fe/H] = -0.08 ± 0.09 (Soderblom k 
Mayor 1993). Soderblom ei al. (1993d) have ob- 
tained ^Li observations for 14 probable members of 
the UMaG which have high quality colors. There 
are few observations of ^Li in cluster stars which are 
older than the Hyades. The cluster NGC 752 has 
an age of 1.7 Gyr (Barry et al. 1987) and there are 
12 ^Li detections and 5 upper limits (Hobbs & Pila- 
chowski 1986a) for this cluster which has solar metal- 
licity. Finally, M67 has an age which is roughly solar, 
4.0^0 5 Gyr (Demarque, Green & Guenther 1992) and 
a similar metallicity to the Sun, [Fe/H] = — 0.04±0.12 
(Hobbs & Thorburn 1991). There are 9 ^Li detections 
and 10 upper limits available for M67 (Hobbs & Pi- 
lachowski 1986b; Spite et al. 1987). We have elected 
not to use the data from the cluster a Per, due to 
difficulties with patchy reddenning and membership 
(Balachandran 1988). 

3. Standard Models 

In order to explore the effect of metallicity and 
overshoot on the models, a number of ^Li isochrones 
were constructed using standard physics (models GD, 
LA and PA from Table 1, Paper I). For each set of pa- 
rameters we are interested in, a series of stellar mod- 
els with masses 0.6, 0.7 ... 1.3 Mq are evolved. The 
basic physics in these models are described in Pa- 
per I. In brief, the energy producing nuclear reaction 
rates are from Bahcall & Pinsonneault (1992); reac- 
tion rates for the ®Li, ^Li and ^Be are from Caugh- 
lan & Fowler (1988); the high temperature opacities 
from Iglesias & Rogers (1991); the low temperature 
opacities (below lO'* K) are from Kurucz (1991); while 
the surface boundary conditions are determined using 
Kurucz model atmospheres (Kurucz 1992; Howard, 
1993). The fit between the stellar model and the at- 
mosphere is made at an optical depth of r = 2/3. 
For temperatures above 10® K, a relativistic degener- 
ate, fully ionized equation of state is use. Below 10® 
K, the single ionization of ^H, the first ionization of 



the metals and both ionizations of '^He are taken into 
account via the Saha equation. 

Theoretical ^Li isochrones are constructed for a 
given age by linearly interpolating the ^Li depletion 
and Teff to the desired age for each individual mass. 
We shall refer to models with different masses but 
similar ages by the two letter code which is used in 
Table 1, Paper I to describe the parameters used in 
the construction of the solar models. Thus, when 
the term 'model JF isochrones' is used, we mean 
the isochrones which were constructed from models 
evolved using the parameters of the JF solar model. 

In order to compare the ^Li destruction isochrones 
to the observations, it is necessary to assume an ini- 
tial ^Li abundance. Observations of ^Li in T-Tauri 
stars are very difficult and produce a wide range 
of abundances, from logN(Li) = 2 ~ 4, where 
logN(Li) = log(^Li /H)-M2 (Basri, Martin & Bertout 
1991). The non-LTE work of Magazzu, Rebolo. and 
Pavlenko (1992) appears to have settled the ques- 
tion of the ^Li abundance in T-Tauri stars; these au- 
thors found logN(Li) = 3.2 ± 0.3 as the abundance 
of T-Tauri stars. The abundance of ^Li in the old- 
est meteorites is logN(Li) = 3.31 ± 0.04 (Anders & 
Grevesse 1988). Observations of ^Li in the interstellar 
medium yield logN(Li) = 3.53 (Lemoine et al. 1993) 
with an likely error of ±0.5 dex (Hobbs, 1993) Thus, 
the initial ^Li abundance of our models is not a well 
determined quantity; it is likely to be in the range 
logN(Li) = 3.3 — 3.6. We have chosen the initial 
^Li abundance which yields the optimal fit to all of 
the cluster data and which lies within the range spec- 
ified above. For models with non-solar metallicities, 
we have scaled the solar heavy element mixture using 
[Fe/H] = log^ZlZo)) and calculated the '^He abun- 
dances by assuming AY = 2.5A2'. 

Figure 1 demonstrates the sensitivity of the models 
to the amount of overshoot below the convection zone. 
Changing the overshoot by 0.05 pressure scale height 
(Hp) has a noticeable effect on the cool models. Mod- 
els with no overshoot, or an overshoot of 0.02 pressure 
scale height, provide a good match to the mean trend 
of the Pleiades observations. However, the models do 
not account for the large spread in ^Li abundances 
which is observed below 5300 K. The standard models 
do not match the mean trend observed in the Hyades 
below 6000 K. The models are under-depleted with 
respect to the majority of the observations. This is 
a clear indication that mixing is active in stars be- 
tween 70 and 700 Myr. Note that the standard models 
do match the observed Lt abundances m the tidally 
locked binaries. As rotational mixing is unlikely to oc- 
cur m tidally locked binaries, this is evidence that the 
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Fig. 1. — The effect of overshoot on the Li destruc- 
tion in standard models is demonstrated by compar- 
ing ^Li destruction isochrones from models PA (Hp = 
0.00), LA (Hp = 0.02) and GD (Hp = 0.05) to the 
observations in the Pleiades and Hyades. The initial 
^Li abundance has been taken to be logN(Li) = 3.3. 
The observations are shown as various symbols, while 
the results of the model calculations are given as lines. 
The Pleiades models have [Fe/H] = 0.0, while the 
Hyades models have [Fe/H] = -|-0.1. 

additional mixing mechanism which operates between 
the ages of the Pleiades and the Hyades is related 
to rotation. The mismatch between standard mod- 
els and observations becomes larger as we go to older 
stars. This is shown in Figure 2 where the theoretical 
isochrones are compared to observations in M67 and 
NGC 752. It is clear that standard models do not de- 
plete enough ^Li to account for the observations. We 
note that a recent paper by Swenson et al. (1994) was 
able to match the observed ^Li abundances in single 
Hyades stars using standard models with enhanced 
oxygen abundances. However, as noted by Swenson 
et al. (1994), all outlying points represent spectro- 
scopic or tidally locked binaries. 

When comparing the theoretical isochrones to the 
observations, it is important to remember that the 
metallicities of the clusters are somewhat uncertain. 
Figure 3 demonstrates the dramatic effect changing 
the metallicity has on the theoretical isochrones. It is 
clear that a change in metallicity of 0.05 dex (which is 
within the 2a errors of the observations) can change 
considerably whether or not we judge a model to pro- 



Fig. 2. — Standard models with different amounts of 
overshoot are compared to Li observations in M67 and 
NGC 752. The initial ^Li abundance has been taken 
to be logN(Li) = 3.3. The observations are shown as 
squares and triangles, and the theoretical isochrones 
are given as lines. The ^Li abundance of the Sun is 
shown by the solar (©) symbol in the M67 graph. The 
NGC 752 isochrones have an age of 1.7 Gyr, while the 
M67 isochrones have an age of 4.0 Gyr. 

vide a good fit to the data. 

It is clear from Figures 1-3 that the amount of 
^Li destruction in stellar evolutionary models is sensi- 
tive to the amount of overshoot included in the mod- 
els and the exact metallicity of the models. Changes 
in metallicity of the magnitude we have studied here 
(A[Fe/H] = 0.25) produce dramatic effects for stars 
with effective temperatures below 6000 K. Chang- 
ing the extant of the overshoot layer by 0.05 pres- 
sure scale-height leads to noticeable effects for mod- 
els with Teff ^ 5500. Thus, it is the cool stars 
which are most affected by the uncertainties in the 
observed metallicity and extent of the overshoot layer. 
Changes in [Fe/H] which are within the 2a errors of 
the observations can turn an unacceptably bad fit 
into a good fit. In addition to the problems with the 
theoretical isochrones, the observations have random 
and possible systematic errors associated with them. 
Due to this extreme sensitivity, we have decided not 
to perform a least squares fit of the ^Li destruction 
isochrones to the data. To do such an analysis would 
require us to put error bars on the ^Li destruction 
analysis, an endeavor which would require explor- 
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Fig. 3. — Standard models with varying metallici- 
ties and an overshoot of 0.05 pressure scale height 
(model GD) are compared to Li observations in the 
Hyades and the Pleiades. The initial ^Li abundance 
has been taken to be logN(Li) = 3.3. The obser- 
vations are shown as various symbols, and the the- 
oretical isochrones are given as lines. The Pleiades 
has [Fe/H] = -0.034 ± 0.024 while the Hyades has 
[Fe/H] = 0.127 ± 0.022. 

ing an enormously large parameter space. Instead 
of a least squares analysis, we have merely chosen 
an initial ^Li abundance (in the range logN(Li) = 
3.3 — 3.6) which yields the optimal 'eye' fit to the 
Pleiades and Hyades ^Li data. When examining the 
^Li destruction isochrones it is important to recall the 
uncertainties in the depletion of ^Li in the cool stars. 
It is clear however, that no change in metallicity or 
overshoot could make the standard models simultane- 
ously fit the observed ^Li abundances in the Pleiades, 
Hyades, NGC 752 and M67. 

4. Pure Microscopic Diffusion Models 

Models which include microscopic diffusion differ 
from standard models in that the diffusion of ^Li out 
of the surface convection zone leads to a depletion 
of ^Li at the surface on the main sequence. The 
models use the microscopic diffusion coefficients from 
Michaud & Proffitt (1993). The coefficients are un- 
certain at the ~ 15%, for this reason, an adjustable 
constant has been introduced which is set to 0.8 
or 1.0. In Figure 4, ^Li observations in the Pleiades 



Fig. 4. — Models which include microscopic diffu- 
sion but do not include rotational mixing, are com- 
pared to Li observations in the Hyades and the NGC 
752. The initial ^Li abundance has been taken to be 
logN(Li) = 3.5. The observations are shown as vari- 
ous symbols, and the theoretical isochrones are given 
as lines. The isochrones shown are from models GD 
(Hp = 0.02), XB (Hp = 0.24) and XC (Hp = 0.24). 

and NGC 752 are compared to models which include 
microscopic diffusion but not rotational mixing. The 
model with an overshoot of 0.02 pressure scale height 
(model ND) provides a reasonable fit to the mean 
trend in the Pleiades, but does not match the ob- 
servations in NGC 752. The ^Li depletion in models 
which include microscopic diffusion but not rotational 
mixing is very similar to the standard models. This 
is due to the long time scale of microscopic diffusion 
in the low mass stars which are studied here. Thus, 
these models have the same problem as the standard 
models do - there is not enough ^Li depletion on the 
main sequence. 

Two models with large amounts of overshoot (0.24 
and 0.25 pressure scale height) are also shown in Fig- 
ure 4. The amount of overshoot in these models pro- 
vides a good fit to the observed ^Li depletion in the 
Sun (models XB and XC). It is clear that such mod- 
els do no fit the observations in the Pleiades; there 
has been far too much pre-main sequence depletion. 
The match to the observations in NGC 752 is rea- 
sonable. This is not too surprising, as we have cali- 
brated the models to match the Sun at an age of 4.55 
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Gyr. Large amounts of overshoot has been suggested 
by Ahrens, Stix & Thorn (1992) as an explanation 
for the ^Li depletion in the Sun. This suggestion is 
only viable if we assume that the Pleiades (and other 
young clusters, such as a Per, whose stars has similar 
^Li abundances to the Pleiades stars) have anoma- 
lously small amounts of overshoot compared to the 
older clusters NGC 752 and M67 and the Sun. This 
is clearly an unacceptable explanation, so we can rule 
out large amounts of overshoot as a reasonable expla- 
nation for the solar ^Li depletion. 

5. Models with Microscopic Diffusion and 
Rotation 

5.1. Parameter Variations 

In this section, the effects of microscopic diffusion 
and rotational mixing are studied by comparing mod- 
els with different parameters to observations of the 
Pleiades and Hyades stars. In order to make the 
inter-comparisons between different models straight- 
forward, we will present isochrones that have the same 
metallicity and age. For the Pleiades, all isochrones 
have an age of 70 Myr and a solar metallicity while 
the Hyades isochrones have an age of 700 Myr and 
[Fe/H] = -FO.IO. 

The details of the rotational mixing prescription 
are found in Paper I, here we provide a brief sum- 
mary. The rotational mixing coefficients are taken 
to be the product of a velocity estimate and the ra- 
dius (Zahn 1993). The velocity estimates are similar 
to those used by Pinsonneault ti al. (1989) and in- 
clude meridonal circulation, dynamical shear instabil- 
ities, the Solberg-Hoiland instability, the Goldreich- 
Schubert-Fricke instability and the secular shear in- 
stability. The rotating models include a general pa- 
rameterization for the loss of angular momentum at 
the surface due to magnetic stellar winds. The for- 
mulation is similar to that given by Kawaler (1988), 
except that a saturation level, Wcrit has been intro- 
duced. If the surface angular velocity is below Wcrit 
then the rate of angular momentum loss is that given 
by Kawaler (1988). Above Wcrit, angular momentum 
loss occurs at a reduced level compared to that given 
by Kawaler (1988). In addition to Wcrit, there are 
four other free parameters in the rotating models: the 
initial rotation velocity given to the calibrated solar 
model (Vi); the power law index for in the angular 
momentum loss law (N); a constant which multiplies 
the diffusion coefficient for the Goldreich-Schubert- 
Fricke instability (fasp), and which characterizes 
the sensitivity of the rotational mixing to gradient in 
the mean molecular weight. 



In order to construct the Li depletion and ro- 
tation velocity isochrones, a series of stellar models 
(masses 0.60, 0.70, ... 1.30 Mq) were calculated with 
[Fe/H] = 0.0 and -|-0.10. For each mass and metallic- 
ity, 3 models were constructed, with initial rotation 
velocities of 10, 30 and 50 km/s. This spans the range 
observed in T-Tauri stars (Bouvier 1991, Bouvier ei 
al. 1993), with the majority of T-Tauri stars hav- 
ing rotation velocities below 20 km/s. For each ini- 
tial rotation velocity, the ^Li destruction and rotation 
rate isochrones were produced in the same manner as 
was done for the standard models. The Radick ei al. 
(1987) rotation period data for the Hyades includes B- 

V colors, not effective temperatures. In order to com- 
pare the models to the observations of rotation rates, 
the effective temperatures have been converted to B- 

V colors using an updated version of the RYI color 
calibration (Green, Demarque & King 1987; Green 
1988). 

Surface rotation and ^Li depletion isochrones have 
been calculated for nearly all of the models presented 
in Table 1, Paper I (Chaboyer 1993). Here, only a 
selected subset of the models will be presented. We 
have selected some models which illustrate the sensi- 
tivity of the models to the various parameters, and 
present a model which best represents the observed 
^Li depletion and rotation velocity patterns observed 
in young cluster stars. We first consider an = 0.10 
model, UU. These models have an overshoot of 0.05 
pressure scale height, N = 2, Wcrit = 3.0 x 10® 
and the microscopic diffusion coefficients are given by 
Michaud & Proffitt (1993). The models are compared 
to observations in Figures 5-6. It is clear from these 
figures that larger initial rotation velocities lead to 
larger ^Li destruction. Figure 5 shows that in model 
UU, the cool Pleiades stars are over-depleted in ^Li, 
while the hot stars are over-depleted in the Hyades. 
The surface rotation isochrones are shown in Figure 6 
with the Pleiades and the Hyades observations. When 
examining Figure 6 (and all subsequent rotation rate 
plots) it is important to realize that the Pleiades ob- 
servations are v sin i data. No attempt has been made 
to correct for inclination effects. The fastest rotators 
have likely been observed with sin i = 1. The aver- 
age statistical factor which converts v sin i data to v 
is 4/7r. The comparison to the rotation periods in the 
Hyades is relatively straightforward, as there are no 
inclination effects to worry about. Model UU fails to 
match the rapidly rotating cool stars in the Pleiades 
and the hot stars in the Hyades. The failure of model 
UU isochrones to fit the observed ^Li abundance in 
the hot stars in the Hyades and to match the rota- 
tion velocities of the cool stars in the Pleiades leads us 
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Fig. 5. — Models UU ^Li destruction isochrones are 
compared to ^Li observations in the Hyades and the 
Pleiades. The initial ^Li abundance has been taken 
to be logN(Li) = 3.5. The observations are shown 
as various symbols. The theoretical isochrones are 
given as lines. The legend identifies the initial rota- 
tion velocity given to the models. The Pleiades mod- 
els have a solar metallicity, while the Hyades models 
have [Fe/H] = -^0.10. 

to reject the UU models as viable stellar models. The 
other models with = 0.1 presented in Table 1, Pa- 
per I did an equally poor (or worse) job of matching 
the observations. 

A number of calibrated models have been cal- 
culated which have with = 0.01, overshoot of 
0.02 pressure scale height, and microscopic diffu- 
sion coefficients multiplied by 0.8. The match be- 
tween the most promising models and the Hyades 
and Pleiades observations are shown in figures 7 - 
11. The JF ^Li isochrones (with /gsf = 10, Wcrit = 
3.0 X 10~® *~^) shown in Figure 7 are less depleted 
than model UU. However, the JF isochrones still do 
not match the hot stars in the Hyades; they are over- 
depleted. JF rotational velocity isochrones are similar 
to the UU cases already shown and yield a poor fit to 
the observations. 

The VN isochrones shown in Figures 8 and 9 have 
/(jgp = 1.0, and Wcrit = 1-5 x 10~® provide a bet- 
ter fit to the data than the previous cases we have 
shown. It is clear that changing the efficiency of rota- 
tional mixing and angular momentum loss has an im- 
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Fig. 6. — Model UU surface rotation isochrones are 
compared to v sin i observations in the Pleiades and 
rotation period observations in the Hyades. The ob- 
servations are shown as various symbols. The theo- 
retical isochrones are given as lines and are identified 
by the initial rotation velocity given to the models. 

portant effect on the models. The VN ^Li isochrones 
show very little dispersion at the age of the Pleiades, 
in contrast to the observations. In addition, the 
fastest rotating models deplete the most Li, in con- 
trast to the observations. However, the match to the 
Hyades is quite good. In particular, the depletion of 
^Li in hot stars (Tgff ^ 5500) matches the observed 
^Li abundances. The surface rotation isochrones are 
shown in Figure 9. Although there are still a few 
very rapidly rotating late type stars in the Pleiades 
which do not fit the models, the overall match to the 
observations is better than the previous cases. The 
predicted Hyades rotational periods match the obser- 
vations well for stars with B — V ^ 0.6. Stars hot- 
ter than this still rotate too fast with respect to the 
models. Overall, it is clear the the VN isochrones rep- 
resent the best match to observations of the models 
considered thus far. Model WE has identical param- 
eters to model VN, except it was assumed that the 
solar ^Li depletion was 100 for model WE (it is a fac- 
tor of 200 in model VN). Although they are not shown 
here, the model WE ^Li isochrones are quite similar 
model VN. This indicates that our models are not too 
sensitive to the exact ^Li depletion assumed for the 
Sun. 
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Fig. 7. — Model JF ^Li destruction isochrones are 
compared to ^Li observations in the Hyades and the 
Pleiades. The initial ^Li abundance has been taken 
to be logN(Li) = 3.5. The observations are shown 
as various symbols. The theoretical isochrones are 
given as lines. The legend identifies the initial rota- 
tion velocity given to the models. The Pleiades mod- 
els have a solar metallicity, while the Hyades models 
have [Fe/H] = -^0.10. 
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Fig. 8. — Model VN ^Li destruction isochrones are 
compared to ^Li observations in the Hyades and the 
Pleiades. The initial ^Li abundance has been taken 
to be logN(Li) = 3.5. The observations are shown 
as various symbols. The theoretical isochrones are 
given as lines. The legend identifies the initial rota- 
tion velocity given to the models. The Pleiades mod- 
els have a solar metallicity, while the Hyades models 
have [Fe/H] = -^0.10. 
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Fig. 9. — Model VN surface rotation isochrones are 
compared to v sin i observations in the Pleiades and 
rotation period observations in the Hyades. The ob- 
servations are shown as various symbols. The theo- 
retical isochrones are given as lines and are identified 
by the initial rotation velocity given to the models. 



3.0 
2.0 



o 1.0 



0.0 

3.0 
2.0 



a> 

o 1.0 



0.0 
-1.0 



Hyades 




Pleiades 



7000 



Te«(K) 



Fig. 10. — Model TJ ^Li destruction isochrones are 
compared to ^Li observations in the Hyades and the 
Pleiades. The initial ^Li abundance has been taken 
to be logN(Li) = 3.5. The observations are shown 
as various symbols. The theoretical isochrones are 
given as lines. The legend identifies the initial rota- 
tion velocity given to the models. The Pleiades mod- 
els have a solar metallicity, while the Hyades models 
have [Fe/H] = -^0.10. 
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Fig. 11. — Model TJ surface rotation isochrones are 
compared to v sin i observations in the Pleiades and 
rotation period observations in the Hyades. The ob- 
servations are shown as various symbols. The theo- 
retical isochrones are given as lines and are identified 
by the initial rotation velocity given to the models. 

Finally, the ^Li depletion model TJ isochrones are 
compared to observations in the Pleiades and the 
Hyades in Figure 10. Model TJ has N = 1.8 and 
/(jgp = 10, otherwise it is identical to model VN. It is 
clear that this parameter set (as do most of the previ- 
ous ones) depletes too much ^Li in hot stars compared 
to the Hyades observations. In addition the rotation 
periods are too long compared to the Hyades observa- 
tions (Figure 11). Hence, model TJ is not compatible 
with the observations. This demonstrates the impor- 
tance the GSF diffusion coefficients have in determin- 
ing the ^Li depletion and surface rotation rates of the 
models. 

5.2. Optimal match to the Observations 

In the previous section, the age and metallicity of 
the isochrones remained fixed as we compared vari- 
ous parameter sets to the observations in the Pleiades 
and the Hyades. It is clear from that study that the 
model VN isochrones were the best fit to the obser- 
vations. An optimal match to the observations in the 
Pleiades, Hyades, Praesepe, UMaG, NGC 752 and 
M67 is made by varying the age and metallicity within 
the 1(7 errors. A comparison to ^Li observations in the 
Praesepe and Hyades is shown in Figure 12. Figure 



o 



3.0 
2.0 
1.0 

0.0 

3.0 

2.0 
1.0 
0.0 



o 



Hyades 




O Single Stars 
X Binaries 

□ Tidally Locl<ed Binaries 
V Upper Limils 




)Singie Stars 
( Binaries 
VUpper Limits 
V, = lOkm/s 



= 30 km/s 
= 50 km/s 



6500 



6000 



5500 



5000 



Fig. 12. — Model VN ^Li destruction isochrones are 
compared to ^Li observations in the Hyades and the 
Praespes. The initial ^Li abundance has been taken 
to be logN(Li) = 3.5. The observations are shown as 
various symbols. 

13 shows the fit to ^Li observations in the Pleiades, 
UMaG, NGC 752 and M67. From these two figures, 
it is clear that our models are capable of reproduc- 
ing in detail the time and mass dependence of the 
observed ^Li abundances. However, the fit to the 
Pleiades does reveal a few problems. The dispersion 
(at a given Tgff ) in the ^Li abundance is much smaller 
in the models than in the observations. In addition, 
there are a few stars (near 5000 K) which lie well 
above the isochrones, indicating that some mecha- 
nism is stopping the ^Li depletion in these stars which 
is not present in the models. The rotational velocity 
isochrones for the best fit cases are very similar to 
those already shown for model VN (Figure 9) and are 
not repeated here. 

6. Summary 

The observed lithium depletion pattern in low mass 
stars is strong evidence for physical processes not in- 
cluded in the standard stellar model. In this paper, we 
have considered the hypothesis that mixing induced 
by rotational instabilities and microscopic diffusion 
is responsible for the observed lithium depletion pat- 
tern. The empirical data requires a mixing mecha- 
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Fig. 13. — Model VN ^Li destruction isochrones are compared to ^Li observations in the Pleiades, UMaG, NGC 
752 and M67. The initial ^Li abundance has been taken to be logN(Li) = 3.5. The observations are shown as 
various symbols. The theoretical isochrones are given as lines. The ages and metallicities have been varied to yield 
an optimal fit to the data. The isochrone ages are given on the plot. The metallicities used are [Fe/H] = —0.05 for 
the Pleiades, and [Fe/H] = 0.00 the others. 



nism which allows for the existence of a dispersion in 
abundance at fixed Tgff in clusters and depletion at 
a declining rate with increased MS age. These are 
natural features of rotational mixing. The observed 
surface rotation rates can be used to constrain the 
angular momentum evolution of low mass stars, espe- 
cially angular momentum loss rates. If the rotational 
mixing paradigm is accepted, then the surface lithium 
abundances can be used to constrain the time scales 
for mixing and internal angular momentum transport. 

The available data on rotation velocities and ^Li 
abundances in young clusters has been examined in 
order to constrain the various free parameters in 
our models. Standard models do not match the 
^Li abundances observed in older clusters such as 
NGC 752 (age 1700 Myr) and M67 (age 4000 Myr). 
However, the standard models do match the mean 
trend in the ^Li abundances observed in the young 
clusters such as the Pleiades (age 70 Myr). This 



is clear evidence that some form of mixing occurs 
in the radiative regions of low mass main sequence 
stars. The addition of microscopic diffusion does not 
solve this problem. In the low mass stars studied here 
(M < 1.3 Mq), the time scale for microscopic diffu- 
sion of ^Li is extremely long. Hence, the ^Li depletion 
in models which include microscopic diffusion are very 
similar to standard models. 

We have constructed a large number of stellar mod- 
els which include microscopic diffusion and rotation- 
ally induced mixing. Due to the uncertainties in the 
time scale for mixing of these instabilities and our 
primitive knowledge of angular momentum loss from 
the surface of a star, there are a number of free param- 
eters associated with the models. When rotational 
mixing is included in the models, the match to the 
^Li abundances in older clusters and the Sun becomes 
quite good. By fine tuning the parameters, it is pos- 
sible to obtain a good fit to the mean ^Li abundances 
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in clusters spanning a wide range in age. 

The ^Li abundance measurements and surface ro- 
tation observations available for the Pleiades and the 
Hyades put tight constraints on the models. These 
observations rule out most of the parameter com- 
binations we have attempted. The best match to 
the observations is obtained using model VN and is 
shown in Figures 9, 13 and 12. The models do a 
good job of matching the mean trend in the observed 
^Li abundances in clusters with ages varying from 70 
Myr to 4000 Myr. However, there are still some prob- 
lems associated with the Pleiades models. In partic- 
ular, there are a number of stars (with Tgff < 5300 
K, M ^ 0.9 Mq, many of which are fast rotators) 
which have higher ^Li abundances than is allowed by 
the models. It is clear that we are not modelling some 
mechanism which is halting the ^Li depletion in cer- 
tain stars. 

The VN isochrones are a good fit to the observed 
rotation periods in the Hyades. The observed t; sin i 
rotation velocities in the Pleiades include a num- 
ber of rapidly rotating cool stars (B — V > 0.9, 
M ^ 0.9 Mq) which our models are unable to re- 
produce. The fact that the stars which are under- 
depleted in ^Li are in the same mass range is intrigu- 
ing. It would appear that certain low mass stars have 
not suffered much angular momentum loss during the 
pre-main sequence phase. This could lead to a smaller 
amount of rotationally induced mixing within these 
stars, which would explain the population of stars 
which have higher than normal ^Li abundances at the 
age of the Pleiades. There are no such outliers ob- 
served in the Hyades, and so these stars must spin 
down between 70 and 700 Myr. 
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